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Abstract 

Motivated by the Minimal Dark Matter scenario, we consider the annihilation into gamma rays of 
candidates in the fermionic 5-plet and scalar 7-plet representations of SU{2)l, taking into account 
both the Sommerfeld effect and the internal bremsstrahlung. Assuming the Einasto profile, we 
show that present measurements of the Galactic Center by the H.E.S.S. instrument exclude the 
5-plet and 7-plet as the dominant form of dark matter for masses between 1 TeV and 20 TeV, 
in particular, the 5-plet mass leading to the observed dark matter density via thermal freeze-out. 
We also discuss prospects for the upcoming Cherenkov Telescope Array, which will be able to 
probe even heavier dark matter masses, including the scenario where the scalar 7-plet is thermally 
produced. 
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I. INTRODUCTION 


While there is solid gravitational evidence for the existence of dark matter (DM), its 
actual nature still eludes us. The most popular hypothesis is that the DM is made up 
of massive and stable thermal relics of the Early Universe [1], a scenario that points to 
DM particles with weak interactions. The existence of such DM candidates is predicted by 
many extensions of the Standard Model (SM) of particle physics, like supersymmetry [2H3]. 
However, in recent years, a mixture of observational and experimental anomalies, together 
with theoretical considerations and the lack of supersymmetric signals at colliders, has lead 
to a complementary, more bottom-up approach. Among the plethora of models that have 
been discussed in the literature, some are characterized by their simplicity, as well as by 
their rich phenomenology. A non-exhaustive list includes singlet scalar PH candidates, the 
inert doublet [HHIQ] and the so-called Minimal Dark Matter 

Minimal Dark Matter, which we take as a motivation for the present work, is a particularly 
predictive framework that extends the SM with just one extra weak multiplet of particles, 
one of which is neutral and is a DM candidate. Assuming that the multiplet interacts 
dominantly through gauge interactions, the only unknown parameter for a given multiplet 
is its mass, which may be hxed if one assumes that the DM was thermally produced in the 
Early Universe. In this paper we will take a more agnostic attitude and we will not make 
any assumption about the DM production mechanism but assume that the mass lies in the 
TeV range. 

Of particular interest are the higher multiplets. In the original proposal it was pointed 
out that a fermionic 5-plet and a scalar 7-plet are the highest multiplets one may consider 
that are stable and are compatible with the Landau pole condition mm- While it was 
pointed out in [15] that the 7-plet may be unstable^, we use these two models as benchmark 
scenarios to analyze gamma-ray signals of DM particles belonging to large SU(2) multiplets, 
with spin 0 and with spin 1/2. Thermal relics from higher multiplets tend to have a very 

^ The stability of the fermionic 5-plet is guaranteed by the fact that it has no dimension-5 effective couplings 
to the SM particles. The scalar 7-plet however has such a dimension-5 coupling, suppressed by a scale 
A, in which three scalar 7-plets in a symmetric triplet state couple to two SM scalar doublets. This term 
leads at one loop to the decay of the DM candidate into SM particles, the lifetime being 0(10“® s) for 
A ~ 10^® GeV and Mdm ~ 1 TeV [15]. To prevent such coupling, an extra discrete symmetry must be 
postulated. 
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large mass, typically a few TeV, lying beyond the reach of the LHC. Constraints from indirect 
searches are thus of greatest signihcance for such DM candidates. 

In the present work we focus on the constraints on DM candidates in large SU{2)l 
representations, concretely the fermionic 5-plet and the scalar 7-plet, that stem from the non¬ 
observation of the gamma-rays which are hypothetically produced in their self-annihilation 
in the Milky Way Center. DM candidates in smaller SU{2)l representations have been 
largely discussed in the recent literature, in particular the fermionic triplet, or wino DM, 
so our analysis is complementary to those of Refs. [TnH22] • The motivation for studying 
higher multiplets, like the 5- and 7-plets, is manifold. First, we want to confront these DM 
candidates to the most recent data from the High Energy Spectroscopic System (H.E.S.S.) 
[23] and to the sensitivity of the future Cerenkov Telescope Array (CTA) [2l], along the 
vein of the forecast analysis made in [25]. Second, as emphasized in UBE], since the DM 
mass is much greater than the electroweak scale, indirect signals from these candidates are 
subject to non-perturbative effects [121 IlHl EZ]- The computation of the non-perturbative 
effects is plagued with numerical instabilities mi, which become even more acute for higher 
multiplets. For our work, and in order to obtain reliable results, we have developed a new 
algorithm which cures the numerical instabilities. Lastly, we also include in our analysis the 
annihilations into a process dubbed internal bremsstrahlung, which generates a 

sharp spectral feature in the gamma-ray energy spectrum [281 - 130] (see also m) Including 
the two-to-three annihilation increases the flux into sharp spectral features and, as we will 
argue, even provides, for certain DM masses, the dominant contribution to the total gamma- 
ray spectrum, due to cancellations among the amplitudes contributing to other processes. 

The rest of our article is organized as follows. In SectionjlTjwe briefly review the properties 
of the fermionic 5-plet and scalar 7-plet DM scenarios which are relevant for our analysis. In 


Section we describe our procedure to evaluate the non-perturbative effects on the anni¬ 
hilation process, and in Section |IV] we derive the cross sections for the relevant annihilation 
channels. Then, in Section|V]we calculate constraints on the model using H.E.S.S data from 


the Galactic Center and, in Section VI, a forecast for the CTA reach. Lastly, we present 


our conclusions in Section VH We also include an Appendix describing our algorithm to 
numerically calculate the Sommerfeld enhancement factors. 
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II. THE 5-PLET AND 7-PLET DM CANDIDATES 


We extend the SM particle content by one extra fermionic or scalar SU{2)i multiplet, 
motivated by the Minimal Dark Matter framework m (see also |32] for a summary). Then, 
the Lagrangian reads 

C = CsM + - M)x (fermion) (1) 


or 

C = Csm + ^{\D^x\^ - M^\x\‘^) (scalar). (2) 

In this work we focus on the fermionic 5-plet and scalar 7-plet as benchmarks for large 
multiplets. Both have zero hypercharge Y = 0, so that the 5-plet is Majorana and the 7-plet 
is real. Then, each of those multiplets can be written in terms of their components as 



^ DM^+ ^ 

f DM2+ ^ 


dm2+ 

DM+ 


DM+ 

DM 

for the 5-plet, y = 

DM 

-DM" 


-DM" 

DM^- j 


DM^- 


^-DM^-y 


for the 7-plet. 


(3) 


Here the relative signs have been introduced in order to make the multiplets isospin self¬ 
conjugate. The notation for the multiplet components is self-explanatory. 

Electroweak symmetry breaking induces a mass splitting between the different particles 
in the multiplet. Both for scalars and fermions, weak isospin breaking effects necessarily 
arise at the one-loop level, through the exchange of Z or IT gauge bosons HH, generating a 
mass splitting between the particles of charge Q and Q' which is given by 

Mq — Mq/ Ri (Q^ — Q''^) A , where A = 02 sin^ ~ MeV. (4) 

(Note that A corresponds to the mass difference between the neutral and the charge ±1 
states.) As a result, the neutral particle is the lightest component of the multiplet. For the 
fermionic 5-plet, the dominant source of isospin breaking is the above-mentioned quantum 
effect, since the isospin breaking effects induced by possibly existing non-renormalizable 
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operators are naturally suppressed. This is not a priori the case for scalar candidates, which 
may have renormalizable couplings to the SM scalar doublet The relevant coupling 

is 

£ D A (x'T‘x) . (5) 

where H is the SM scalar doublet and T“ are the generators of the SU{2)l algebra. Nev¬ 
ertheless, {x^T°‘x) is identically vanishing for real scalars P!. which is the case we are 
considering here, therefore the splitting for the scalar 7-plet is naturally dominated by loop 
corrections, as in the fermionic case. 

If the fermionic 5-plet or scalar 7-plet candidate is the dominant form of DM, then 
the only unknown parameter for a given multiplet is the DM mass, which may be hxed 
by requiring that the DM is generated via thermal freeze-out. Taking into account the 
Sommerfeld enhancement in the Early Universe, this requirement leads to a 5-plet mass of 
about 10 TeV and a 7-plet mass of about 25 TeV na. Nevertheless, other scenarios are 
also possible. For instance, for masses below the above mentioned values, thermal freeze- 
out of DM particles only produces a fraction of the total DM density, therefore other DM 
species would be necessary to explain the astrophysical and cosmological data (see, e.g. [S3])- 
Yet, annihilations of DM particles in the Galactic Center could generate observable signals. 
On the other hand, for masses above those values, thermal freeze-out overproduces DM 
particles. However, some dilution mechanism could be at work reducing the DM density 
(see for instance [3Tj for a discussion of some possible alternative non-thermal histories). 
Besides, DM particles might have never reached thermal equilibrium, as it is the case in 
cosmological scenarios with a reheating temperature below the DM mass. In this case, some 
non-thermal mechanism should be postulated to generate the observed DM abundance (or 
a fraction of it). Given our ignorance of the dark sector and of the thermal history of the 
Universe, we will leave both the DM mass and the DM density today as free parameters, 
and we will concentrate on the indirect detection of this class of scenarios. In particular we 
will concentrate on the constraints on DM annihilation into gamma-rays from observations 
of the Milky Way Genter, most notably from searches for sharp spectral features. Since 
the annihilation cross sections are fixed for a given DM mass, the non-observation of a 
gamma-ray excess will lead to constraints on the density fraction of the given DM candidate, 
regardless of any hypothesis about the DM production mechanism. 

In the next section we will discuss in more detail our procedure to calculate the non- 
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perturbative effects in the annihilation in the Milky Way Center of the fermionic 5-plet or 
the scalar 7-plet and which play a central role in correctly assessing the prospects for indirect 
detection of these DM candidates. 

III. THE SOMMERFELD EFFECT 

The fermionic 5-plet and the scalar 7-plet DM particles annihilate via the exchange of 
gauge bosons. Since these are very light compared to the annihilating particles, the long 
range interaction can signihcantly distort the wave function of the non-relativistic particles 
involved in the annihilation. As a result, the perturbative calculation of the annihilation 



5-plet 

7-plet 

V{r) 

'^4(7 -2B 0 \ 

-2B C -3\/2R 

0 -3V2B 0 y 

^ 9C -3 B 0 0 \ 

-3B 4C -5B 0 

0 -5R C -6V2B 

y 0 0 -GV2B 0 y 

D— 9^ „-mwr /^_oA 

“ 47rr^ ^ ~ A-wr A-nr ^ 

y^w+w- 

(2 5 3a/2) 

(3 8 11 6\/2) 


V^e2(22 1 0) 

V^e2(32 22 1 0) 

y^zz 

tan^ W-j/'y 

y^^z 

w 

tan6»iy 


Table I: Potential matrices for the 5-plet and 7-plet representations, as well as row vectors W 
entering the calculation of the annihilation cross sections (c/. Eq. Q). The DM pair state vector 
is (DM2+DM2-, DM+DM-, DM DM) for the 5-plet and (DM^+DM^", DM^+DM^", DM+DM”, 
DM DM) for the 7-plet. 
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cross section breaks down and instead a non-perturbative approach ought to be pursued, 
leading generically to a larger value of the cross section. This is the so-called Sommerfeld 
enhancement, hrst described in Ref. [35] for electron-nucleon scattering at low relative ve¬ 
locities. In this paper, we calculate this effect following closely the formalism presented 
in Ref. [26] and [12], consisting in hrst calculating the non-relativistic potential associated 
to the exchange of gauge bosons and then solving the corresponding Schrodinger equation 
for the two-state wave-functions involved in the annihilation. From the solution one ob¬ 
tains a set of non-perturbative enhancement factors, which are hnally used to calculate the 
annihilation cross section in combination with the perturbative amplitudes. 

For the 5-plet case, due to the exchange of gauge bosons, the dark matter pair states 
DM DM mix with the states DM"*" DM“ and DM^’'“DM^“. Consequently, the matrix g{r) 
introduced in Ref. [26], and which encodes the Sommerfeld effect, is 3 x 3. For the 7-plet 
case, DM pairs additionally mix with the state DM^’'"DM^“, and the resulting g{r) matrix 
is 4 X 4. These matrices satisfy the Schrodinger equation [35] 

^Mv‘^1 — V{r)^ g{r) = 0 , (6) 

where v is the relative velocity of the DM particles, concretely v = 2 x 10“^ in our analysis, 
and V (r) is a central potential, which includes the exchange of gauge bosons among pairs 
of particles and the mass splitting among them. We will consider the s-wave piece of the 
annihilation cross sections both for scalar and Majorana particles, therefore the potential 
matrices V (r) for each SU{2)l representation are the same for scalars and for fermions; their 
explicit form is shown in Table 

Since the DM relative kinetic energy in the Milky Way Center is smaller than the mass 
splitting between the DM and the charged particles of the multiplet, only DM pairs can exist 
at large distances in the scattering process. The boundary conditions for the Schrodinger 
equation Eq. (|^ then read 



^/(O) = 1, 

g(^r) -)■ e(*^"’/2)A/i-4y(oo)/(Mt.2)^ ^ 


^irMvl2 


^ 0 
0 


0 ^ 
0 


when r —)• oo . 


(7) 


V" 


(i_i— 
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Here D is a constant matrix and d = (■ ■ ■ _doo) is the only non-zero row of the matrix 

g{r) after factorizing ont the ont-going wave describing the DM pairs at inhnity. Notice that 
when gange interactions are tnrned-off, d is just a unit vector. Consequently, its components 
can be interpreted as non-perturbative enhancement factors due to the exchange of gauge 
bosons. These components, furthermore, are all real up to a global phase.^ 

We describe in Appendix|^a novel procedure for solving Eq. (|^ which cures the numerical 
instabilities which plague the calculations when A/M -C 1 |T^. Using this method, we 
calculate the Sommerfeld enhancement factors for the 5-plet (7-plet) for masses between 1 
TeV and 30 TeV (75 TeV), which approximately corresponds to three times the mass required 
by thermal production. The resulting values are shown in Figs. [^and[^ The solid (dashed 
line) indicates the values of the mass where the enhancement factors are positive (negative), 
while the gray band indicates the value of the mass reported in Ref. na that leads to the 
observed DM abundance via thermal freeze-out, and which approximately corresponds to 
M ^ 10 TeV for the 5-plet and M ^ 25 TeV for the 7-plet. As apparent from the plot, and 
as expected for the non-perturbative exchange of light particles, we hnd multiple resonances 
across the TeV scale in both cases. 


IV. ANNIHILATION CROSS SECTIONS 


We consider now the cross sections that are relevant for calculating the gamma-ray spec¬ 
trum from DM annihilations in the Galactic Center. The prompt gamma-ray flux receives 
three main contributions: i) a soft featureless part consisting of secondary photons from the 
decay and fragmentation of the W and Z pairs that are produced in DM annihilations, ii) 
the mono-energetic photons produced in the processes DM DM —)■ 77 and DM DM —)■ jZ, 
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This result follows from the fact that the potential matrices V{r) are not only Hermitian but also real. 
Then, using Eq.Q, it is possible to prove that the following matrices are conserved 

Ji = 9 {rV 9 '(r) - g'{r)'^ g{r), J 2 = i {gir^ g'{r) - g'{r)^ g{r)) . 


Noticing that the boundary conditions Eq. 0 at infinity imply Ji = 0 and J 2 = Mvdld, while at r = 0 
A = g'{0) - ff'(O)^ and J 2 = i{g'{0) - g'iO^), one obtains that the matrix J 2 = Mvdld must be real. 
Finally, since the row vector d can in principle be complex, one concludes that the Sommerfeld factors 
must share a common global phase, which is irrelevant for the evaluation of the annihilation cross sections, 
as follows from Eqs. § or 




doo 



d+- 



di+i- 



Figure 1: Sommerfeld enhancement factors for the fermionic 5-plet. The solid (dashed) line cor¬ 
responds to positive (negative) values of the enhancement factors. The gray band indicates the 
value of the mass reported in Ref. |12) which leads to the observed DM abundance via thermal 
freeze-out. 

and in) the photons prodnced in the internal bremsstrahlnng process DM DM —)■ 

The last two contribntions are particnlarly relevant for indirect dark matter detection since 
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Figure 2: Same as Fig. but for the scalar 7-plet. 


they produce sharp spectral features in the gamma-ray spectrum for energies close to the 
DM mass, which might stand out over the featureless astrophysical background, thus offering 
a powerful probe of dark matter models |2H1 1361 - HT] . 
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The s-wave cross section for DM annihilation into a two-body final state reads 


((jn)v'iV2 ~ 2d* rviV2 


with 




( 8 ) 


where the vector of Sommerfeld enhancement factors d was calculated in Section m and 
Wvi y. vectors given in Table which follow from the tree-level annihilation ampli¬ 

tude of the different pairs into a given gauge boson pair Vi V 2 . Here, ^ = 1 for Majorana 
DM whereas ^ = 2 for scalar DM. 

We show in Fig. |^the DM annihilation cross sections into W^W~ and 77 as a function 
of the DM mass for the fermionic 5-plet and the scalar 7-plet; the cross sections into ZZ 
and 'jZ can be readily obtained from {<Jv)ry^ from the relations 


{o-v)zz 


(^^)77 

tan"^ 6w ’ 


{av)^z 


tan^ 9w 


(9) 


The resonance structure at the TeV scale, originated by the non-perturbative effects, is 
apparent in Fig.|^ Furthermore, we find dips in the two-body annihilation channels, namely 
masses for which the cross sections vanish, and which result from the destructive interference 
between Sommerfeld factors (this behavior, known as Ramsauer-Townsend effect, was also 
noted in [I 6 ]). For instance, the 5-plet cross section for the 77 channel reads 

(H 77 = -§r M+- + 4d2+2-|^ • (10) 


Therefore, for the mass ranges where the Sommerfeld factors have opposite signs, the cancel¬ 
lation between the two terms leads to a reduction of the cross section. An analogous effect 
also occurs for other annihilation channels. These cancellations are particularly signihcant 
for the scalar 7-plet, which presents dips in the W~^W~ and 77 cross sections for values of 
the DM mass comparable to the thermal value; we have checked numerically that the posi¬ 
tion of the dips does not change significantly when the relative velocity of the DM particles 
is increased or reduced. On the other hand, for the fermionic 5-plet we hnd no cancellations 
in the cross sections in the mass range of relevance for thermal production. 

In addition to the two-to-two annihilation channels, in this class of models two-to-three 
annihilations can play an important role in indirect searches. Indeed, and due to the mass 
degeneracy among the components of the DM multiplet, the internal bremsstrahlung gener- 
ically leads to a bump in the gamma-ray energy spectrum close to its kinematical endpoint. 
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Fermionic 5-plet 



M (TeV) 
Scalar 7-plet 



Figure 3: DM annihilation cross section into W~^W (blue line) and 77 (green line) for the fermionic 
5-plet (upper panel) and the scalar 7-plet (lower panel). 


which might stand out over the astrophysical background. The Sommerfeld enhanced inter- 
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nal bremsstrahlung cross section reads 


dE^ 


1287r3M2 



1 

71 


dooM (DMDM ^ W+W--f) 
2 


+ J2di+i-M (DM*+DM*- ^ W+W-j) 

i 


dEy[r+ . 


( 11 ) 


where E^+, are the minimum and maximum energies kinematically accessible to the W 
boson produced in the annihilation, while Ai stands for the tree-level annihilation amplitude, 
which we obtain using FeynArts with an implementation of the fermionic 5-plet and 
scalar 7-plet DM models made with FeynRules |13]. 

In Fig. 1^ we show, for illustration, the main components of the prompt gamma-ray 
spectrum produced in the annihilations of a fermionic 5-plet with M = 10 TeV (left plot) 
and scalar 7-plet with M = 25 TeV (right plot), assuming an energy resolution of 10%. 

We also remark that for small DM relative velocity, as is the case of annihilations at 
the Galactic Center, the Sommerfeld enhancement factors are identical for fermionic and 
for scalar DM, since the potential matrix entering the Schrondiger equation for the relevant 
two-body states, namely those with total charge Q = 0 and total spin S' = 0, is the same for 
a given SU{2)l representation. Therefore, for a given DM multiplet, the annihilation cross 
section for scalar DM is simply a factor of two larger than for Majorana DM, as a result of 
the different value of the parameter ^ in Eqs. (§. (0. The differential annihilation cross 
section for the scalar 5-plet (fermionic 7-plet) DM candidate, as well as their experimental 
limits, can then be straightforwardly derived from the analysis presented in this paper for 
the fermionic 5-plet (scalar 7-plet). 
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V. CURRENT LIMITS FROM H.E.S.S. 


The H.E.S.S. instrument has measured the gamma-ray flux in the direction of the Galactic 
Center in an energy range from 300 GeV to 30 TeV, finding an excellent agreement with the 
background expectations, thus allowing to set constraints on the rate of DM annihilations. 
The expected gamma-ray flux from annihilations in the solid angle AD centered in given 
direction of the sky is given by 


d(j)^ 1 d{av] 

dE^ dE^ 


J. 


( 12 ) 


where Jann is the astrophysical J-factor, calculated by integrating the square of the DM 
density profile of the Milky Way p(r) over the line of sight and the region of interest AD 




/ dD / dsp 
'Afi J\os 


(13) 


In our analysis we employ the gamma-ray flux measured with the H.E.S.S. instrument in 
a target region consisting of a circle of 1° radius centered in the Milky Way Center, excluding 
the Galactic Plane by requiring |6| > 0.3° |2Sl SI]- In order to bracket the astrophysical 
uncertainties in the calculation of the flux, we consider two different DM halo profiles: the 
cuspy Einasto profile 


PEin = P^exp 




1 

[ a 

ArJ 

i 


(14) 


with scale radius = 20kpc and shape parameter a = 0.17, and the cored isothermal 
profile 

Piso = , , sp , (15) 


1 + (r/r,)2 ’ 

where = 5kpc m- Here, ps is chosen such that the DM density at the position of 
the Sun, which we consider located at a distance d© = 8.5 kpc from the Galactic Center, 
reproduces the local DM density po = 0.39GeV/cm^ |18]. Then, the J-factors are Jeiii = 
4.43 X lO^^GeV^ sr/cm^ [H] for the Einasto profile and Ji^o = 3.23 x lO^^GeV^ sr/cm^ for 
the isothermal profile. 

On the other hand, the differential cross section to gamma-rays d{av)^/dE^ consists 
of three parts, as mentioned in the previous section. First, a continuum of gamma-rays, 
which is produced in DM annihilations into W and Z bosons, and which is calculated 
from the product of the cross section into W^W~ or ZZ and the differential number of 
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Figure 4: Differential annihilation cross sections into 77 (magenta line), 7 Z (pink line), W~^W~^ 
(blue line) and continuum from ZZ and W~^W~ (green line), as well as the total differential cross 
section (black line), for a fermionic 5-plet DM candidate with mass 10 TeV (left panel) and a scalar 
7-plet with mass 25 TeV (right panel). We assume an energy resolution of 10% which is typical 
for current gamma-ray telescopes. 


photons produced per annihilation, which we take as dN^/dx = 0.73 x~^'^ , with x = 

E... 1 /M [19]. Second, the mono-energetic photons produced in annihilations into 7 Z and 
77 , which would be observed with H.E.S.S. as Gaussian distributions, normalized to 1 
and 2 respectively, with a width given by the energy resolution of the instrument, which 
is approximately 17% at 500 GeV improving to 11% at ITeV [23] • Finally, the line-like 
contribution from internal bremsstrahlung, given by Eq. (0 , which we also smear out with 
the energy resolution of the instrument. In Fig. |^we show, for illustration, the differential 
cross sections for red a fermionic 5-plet with M = 10 TeV and a scalar 7-plet with M = 
25 TeV, which approximately correspond to the thermal masses of these two DM candidates. 

We now calculate constraints on the fermionic 5-plet and scalar 7-plet DM models which 
follow from the non-observation by H.E.S.S. of a sharp gamma-ray spectral feature nor an 
exotic featureless contribution to the gamma-ray flux. To calculate limits on the DM annihi¬ 
lation cross section into sharp spectral features, we adopt the phenomenological background 
model proposed by the H.E.S.S. collaboration |23|, which is described by 7 parameters. We 
then include the sharp spectral feature from DM annihilation and calculate the 95% G.L. 
one-sided limit on the cross section. On the other hand, we calculate constraints on the 
featureless component from W^W~ 01 ZZ annihilations following [H], which compares the 
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gamma-ray fluxes measured with the H.E.S.S. instrument in a “search region” and in a 
“background region”. The inferred residual flux is consistent with zero, thus allowing to 
derive upper limits on the flux from annihilations. We note that this approach is only valid 
for cuspy prohles, since for cored prohles the signal gets completely subtracted, therefore in 
this case we will only present limits for the Einasto prohle. 

We show the impact of the H.E.S.S. measurements on the fermionic 5-plet and scalar 
7-plet DM models in Figs. and respectively, for the Einasto and isothermal prohles, and 
separating the limits from the non-observation of sharp spectral features (blue solid lines) 
or from an excess in the continuum (green solid lines). We present the constraints in terms 
of the DM fraction, which is given by the square root of the signal normalization factor. 
We highlight the role of the internal bremsstrahlung in probing the fermionic 5-plet and 
scalar 7-plet models by showing also the limit obtained when only the gamma-ray lines from 
annihilations into 77 and 7 Z are included (dotted blue line). As apparent from the plots, 
the limits from lines exhibit peaks, which correspond to the dips in the annihilation cross 
section shown in Fig. and where cancellations among amplitudes occur. On the other 
hand, the annihilations into W^W~'^ are not affected by these cancellations, thus allowing 
to probe DM masses which are otherwise unconstrained by line searches. Furthermore, after 
including the W^W~'y channel, we hnd that for DM masses above ~ 1.2 TeV the limits 
from sharp spectral features are more stringent than those from the continuum, concretely 
by a factor of 1.5 ~ 8.5 for the fermionic 5-plet and 1.5 ~ 11 for the scalar 7-plet. 

Our analysis excludes DM consisting exclusively of fermionic 5-plet or scalar 7-plet neutral 
states in the whole mass range between M = 1 TeV up to M = 20 TeV, assuming the 
Einasto prohle. More concretely, for low masses the limits are very stringent and exclude 
a DM fraction of 7 x 10“'^ (1.5 x 10“"^) for the 5-plet (7-plet) at M = 1.6TeV (2TeV). In 
the case of the isothermal prohle, we hnd some mass ranges where the 5-plet can be the 
dominant component of DM, notably the regions around the 5-plet thermal mass. The 7- 
plet, on the other hand, is excluded as the dominant component of DM for masses between 
1 and 20 TeV, with the exception of a narrow mass window around 12 TeV. Thermal masses 
for the 7-plet are, unfortunately, not accessible to the H.E.S.S. instrument, although they 
may be probed by CTA, as we will discuss in the next section. 
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Fermionic 5-plet, Einasto profile 



M (TeV) 

Fermionic 5-plet, isothermal profile 



M (TeV) 

Figure 5: 95% C.L. limits on the dark matter fraction for the fermionic 5-plet from the non¬ 
observation by H.E.S.S. of sharp gamma ray spectral features, including (solid blue line) and 
neglecting (dotted blue line) the internal bremmstrahlung contribution, as well as from the non¬ 
observation of the continuum gamma-rays from annihilations into W~^W~^ and ZZ (green line), 
assuming the Einasto profile (upper panel) and the isothermal profile (lower panel). 
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Scalar 7-plet, Einasto profile 



M (TeV) 

Scalar 7-plet, isothermal profile 



Figure 6: Same as Fig. but for the scalar 7-plet. 


VI. PROSPECTS WITH CTA 

We consider in this section the sensitivity to the fermionic 5-plet or scalar 7-plet DM 
candidates of the npcoming Cherenkov Telescope Array, an instrnment that will be able to 
measnre the gamma-ray flnx np to energies ~ 100 TeV. In onr analysis we will employ the 
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instrument properties reported in [50] for the array I, a balanced array that can measure 
gamma-rays in a wide energy range from a few tens of GeV up to over 100 TeV. Its effective 
area is larger than 10® m^ above 1 TeV with a resolution better than 10% in this energy 


3 


range. 


Since in these scenarios the limits from continuum gamma-rays are in general worse than 
the limits from sharp spectral features, we will only consider the latter for our CTA forecast. 
In our analysis we closely follow the approach described in |25]. First, we generate mock 
data. The expected number of counts in the energy bin i is given by 



( 16 ) 


where At is the observation time in the region of interest, the width of the ith energy 
bin and <htot the total background in the region of interest. The energy resolution R{E, E') 
and effective area Aefr of the instrument are taken from [50]. In this analysis we use 200 
energy bins per decade in order to resolve the sharp spectral features. Then, from the mean 
number of counts in each energy bin, we generate 200 sets of mock data using Poisson 
distributed random numbers. 

In order to compare the CTA reach with the H.E.S.S. constraints, we consider the same 
region of interest around the Galactic Center and 112h of observation time. For the back¬ 
ground we make use of the parametrization presented in [25], which takes into account not 
only the diffuse gamma-rays and the H.E.S.S. point source in the Galactic Center, but also 
cosmic ray backgrounds. In fact, protons are so abundant in cosmic rays that they constitute 
an important contribution, especially at higher energies, even with good proton rejection 
efficiencies of 10“^ — 0.2. Electrons and positrons produce air showers that are indistinguish¬ 
able from gamma-ray initiated showers, hence they comprise an irreducible background for 
Cherenkov telescopes. Eor our region of interest, electrons are the dominant background at 
low energies. 

The limits are calculated using the sliding energy window technique. This method is 
based on the assumption that the background is locally well described by a power law. The 

^ The CTA collaboration has recently released an updated performance for a different array configura¬ 
tion m- We have checked that, with respect to array I, the corresponding limits on the DM fraction for 
gamma-ray monochromatic lines improve up to a factor of 2.8 for low masses and up to a factor of 1.3 
above 4 TeV. 
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energy range for which this is the case determines the size of the energy window. In our 
analysis we use the largest possible window size of e = 2 |2S], which is suitable for the CTA 
background. Finally, inside these energy windows, we determine the 95% C.L. one-sided 
upper limits using a prohle likelihood analysis. 

The resulting limits for both prohles under consideration are shown in Figs, [^and [^for 
the fermionic 5-plet and scalar 7-plet, respectively. Again, we present the limits in terms of 
the DM fraction, which is obtained from logarithmically averaging the signal normalization 
factor in the 200 mock data realizations, and then taking its square root. We hnd that CTA 
can improve the H.E.S.S. limits on the DM fraction by a factor of 1.2 ~ 3 and close in on 
scenarios with heavier masses. More specihcally, even for the isothermal prohle, CTA will 
be able to exclude 5-plet masses up to 13 TeV and 7-plet masses up to 40 TeV. In particular, 
CTA will be able to exclude the possibility that the 7-plet is produced via thermal freeze- 
out. On the other hand, for the Einasto prohle, DM entirely consisting of 5-plet and 7-plet 
states can be excluded up to 30 TeV and 75 TeV, respectively, which roughly corresponds 
to three times the values of the mass favored by thermal production. 
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Fermionic 5-plet, Einasto profile 



M (TeV) 

Fermionic 5-plet, isothermal profile 



Figure 7: 95% C.L. projected limits on the dark matter fraction for the fermionic 5-plet from 
the non-observation by the upcoming CTA of a sharp spectral feature in the Galactic Center, 
including (solid magenta line) and neglecting (dotted magenta line) the internal bremmstrahlung 
contribution, assuming the Einasto profile (upper panel) and the isothermal profile (lower panel). 
We also show for comparison the current H.E.S.S. limits (solid blue line). 


21 














Dark Matter fraction Dark Matter fraction 


Scalar 7-plet, Einasto profile 



M (TeV) 

Scalar 7-plet, isothermal profile 



Figure 8: Same as Fig. but for the scalar 7-plet. 
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VII. CONCLUSIONS 


In this work we have studied the gamma-ray spectrum generated in the annihilation of 
fermionic 5-plet or scalar 7-plet dark matter particles, motivated by the Minimal Dark Mat¬ 
ter framework. We have remained agnostic about the dark matter production mechanism in 
the Early Universe, hence the dark matter mass is a free parameter in our analysis, which we 
take between 1 TeV and 30 TeV (75 TeV) for the 5-plet (7-plet), which roughly corresponds 
to three times the mass expected for a thermal relic. In this mass range, non-perturbative 
effects are expected to be signihcant and ought to be taken into account. In this work, we 
have also developed a new algorithm to calculate the Sommerfeld enhancement factors in the 
limit of low relative velocities, which cures the numerical instabilities reported in previous 
approaches. 

We have included in our analysis the internal bremsstrahlung, which produces a sharp 
spectral feature in the gamma-ray spectrum as a natural consequence of the small mass 
splitting between the charged and neutral states of the dark matter multiplet. This compo¬ 
nent leads to new opportunities to detect sharp gamma-ray spectral features in this class of 
scenarios, since the cross sections into 77 and 7 Z vanish for certain values of the DM mass, 
due to cancellations among the amplitudes contributing to these channels. 

We have then confronted the fermionic 5-plet and scalar 7-plet dark matter scenarios 
with the most recent Galactic Center observations of the H.E.S.S. instrument, considering 
the continuum of gamma-rays from annihilations into W~^W~ and ZZ, as well as the sharp 
spectral features from monochromatic photons and internal bremsstrahlung. Assuming the 
Einasto prohle, the H.E.S.S. measurements already exclude DM consisting exclusively of 
fermionic 5-plet or scalar 7-plet neutral states in the mass range from 1 TeV up to 20 TeV, 
which includes in particular the mass expected for a thermally produced 5-plet DM. We 
have also discussed the sensitivity of the upcoming Cerenkov Telescope Array to this class 
of scenarios, in particular, the prospects to detect sharp gamma-ray spectral features. We 
hud that CTA can improve the H.E.S.S. reach to the DM fraction in this class of models by 
a factor of 1.2 ~ 3. For the Einasto prohle, dark matter entirely consisting of 5-plet (7-plet) 
states can be probed in the mass range between 1 TeV and 30 TeV (75 TeV), and will be able 
to test these two scenarios assuming thermal DM production. For the isothermal prohle, the 
prospects are poorer but can nonetheless probe large regions of the parameter space where 
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either the 5-plet or the 7-plet constitutes the dominant component of dark matter. 

We conclude highlighting the role of Cherenkov telescopes in probing scenarios with 
heavy dark matter particles, as is the case of scenarios motivated by the Minimal Dark 
Matter framework, and which are challenging to probe with direct detection experiments or 
with collider experiments. 


Appendix A: Numerical Algorithm for Calculating the Sommerfeld Effect Today 


The Sommerfeld enhancement factors are obtained solving Eq. ([^ with the boundary 
conditions given in Eq. Q. Although this is in principle straightforward, the numerical 
calculation of the solution is challenging, specially when A/M ^ 1, due to the particular 
form of the boundary conditions im. Here we introduce a novel algorithm that cures the 
numerical instabilities that plague the calculation of the enhancement factors. 

On the one hand, let us notice that Eq. @ can be rewritten as 

h'{r) + h{rY + M — H(r)^ = 0 , where h{r) = g'{r)g{r)~^. (Al) 

with a boundary condition at large distances given by 

h{oo) = (iMn/2)y/l - W{oo)/ (A2) 


In this expression the eigenvalues are either purely imaginary or negative, so that the solution 
at inhnity is either an out-going or a damping wave, respectively. 

On the other hand, conservation of probability on Eq. 0 can be written as 

{g{r)^g'{f) - g{r)'^g{r)) = 0 , (A3) 

which along Eq. ([^ can be used to prove that 

(A4) 


where d = {■ ■ ■ doo)- 

Therefore, in order to calculate the enhancement factors, instead of solving the second- 
order differential Eq. 0 with boundaries at two different points, it suffices to solve the 
hrst-order differential equation Eq. (Al) with one boundary at inhnity. With this solution. 


we calculate the right-hand-side of Eq. (A4), which must have rank one when the relative 
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velocity of the DM particles is small (more concretely, when the DM relative kinetic energy 
is smaller than the mass splitting among the particles in the multiplet). Finally, we calculate 
the corresponding eigenvector to determine the enhancement factors (up to a global phase). 
Notice that this method cannot be used in the Early Universe, since pairs of charged states 
can exist at inhnity at hnite temperature, thus yielding a matrix d^d with rank higher than 
one. 

Note Added 

During the preparation of this work, we became aware of another group investigating 
gamma-ray signals of Minimal Dark Matter scenarios [52], and more recently, about the 
analysis carried out in Ref. [S3], which also addresses the indirect detection of heavy dark 
matter multiplets. 
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